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ABSTRACT 

Context. It is known that the Canary Islands are normally affected by dominant winds flowing from north-northeast, that in some 
meteorological conditions, can transport sand from the Sahara desert at high altitude. The dust may affect the efficiency of the 
telescopes and decrease the transparency of the sky. 
Aims. To maximize the scientific return of the telescopes located at the Observatorio del Roque de Los Muchachos (ORM), we present 
an analysis of the atmospheric dust content and its effects on astronomical observations. Than B, V and / dust aerosol astronomical 
extinction are derived. 

, Methods. Using a 5-year series database of data taken from a dust monitor located inside the Telescopio Nazionale Galileo (TNG) 

dome, we computed mean hourly and daily values of the dust content as measured with a four-channel dust monitor. 
■ Results. We detected particles of 0.3, 0.5, 1.0, and 5.0 fim. Furthermore, using a power law we derived the content of 10.0 fim 

' particles. We found a typical local dust concentration ranging from 3 • 10 6 particles per cubic metre at 0.3 /jm, to 10 3 at 5.0 fim 



and 10 at 10.0 /im, increasing up to 3 orders of magnitudes during the dust storms, with a relatively higher increase of 1.0, 5.0, and 
10.0 fim particles. The number of local dust storm events is the same in the local winter and summer, but the average background 
and storm-related increases in the dust concentration in summer are significantly higher than in winter. In a uniform approximation, 
during the dust storms, an average height of the dust layer of 2.5 km above the telescope is inferred. 

Conclusions. During the sand storms, La Palma Island is affected by an almost uniform layer extending up to 5 km above the sea 
level. The visible extinction is dominated by particles at 1.0, 5.0 and 10.0 fim. In agreement with the results from Carlsberg Automatic 
■ Meridian Circle (CAMC), we find a typical extinction of about 0.2 mag airmass -1 during dust storms. 

Key words, site testing, atmospheric extinction, dust pollution 

X 

' 1. Introduction a knowledge of the amount of the dust and, if possible, of the 

resulting extinction is very important for optimizing the per- 

The Telescopio Nazionale Galileo (TNG) is one of the interna- f ormance f both instruments and telescopes and for obtaining 

tional telescopes located at the Observatorio del Roque de Los re liable photometric measurements. Moreover, the knowledge 

Muchachos (ORM) on La Palma Island (Canary Arcipelago). of the grain sizes and their distribution at the mirror level is 

Due to the proximity of the North African continent, the winds fundamental because the spectral index of the extinction depend 

coming from the central side of the Sahara desert carry a large on the size of grains In fact> grains having dimensions larger 

amount of dust in suspension that can reach high altitudes and ^ me ad opted wavelength produce a neutral extinction, while 

increase the atmospheric dust concentration. sub-micron particles produce an extinction in the same band, 

Whittet et al. (U987|l pointed that this wind action removes with resp ect to the adopted wavelength. 

approximately 2 • 10 8 tonnes of dust every year from Sahara The first detailed study of the impact ofCali ma on astronomical 

desert, and they can be carried out to altitudes up to 6 km. This observations in La Palma is in Murdin ([19851 who found 

seasonal action can reach the south of the Canary Islands, and a seaS onal trend in atmospheric transmission, ranging from 

extended dust clouds often appear above the sky of La Palma, 90% of nights classified as good (i.e., nearly aerosol-free), 

originating a phenomenon locally named la Calima. between October and March, to 58% during July, August, 

This natural event is often discuss concerning the quality and September. Another interesting result from Murdin is the 

of astronomical observations obtained at ORM such when atm0 spheric extinction measured in a site-test campaign in La 

quantify the contribution of astronomical extinction. Therefore Palma A Murdin summarisation of studies from other authors 

shows that 78% of the clear nights in the months starting from 

Send offprint requests to: G. Lombardi February to September 1975 and 83% from December 1974 to 
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November 1975 have extinction less than 0.3 mag airmass . 
To improve the TNG meteorological station, the TNG site 
group put in operation a dust monitor. The first preliminary 
results can be found in Porceddu et al. 120021 and Ghedina et 
al. 120041 The primary motivation behind this measurement 
is to have a continuous characterisation of the air conditions 
above the primary mirror for several reasons. The first one is 
to compare dust measurements and atmospheric extinction to 
know which level of dust may affect astronomical observations. 
The second is to have better maintenance. In fact we know that 
the deposition of the dust over the optical surfaces reduces both 
the optical throughput (laboratory measurements of the primary 
mirror reflectivity show a decrease from 99% reached in the 
case of fresh aluminum to 70% in the case of a dusty surface) 
and the signal-to-noise ratio of the images. The third one is 
for safety reasons and for coating lifetimes, when humidity 
condensation sticks the dust on the mirror. The knowledge of 
the level of dust counts helps us to prevent damage to the optical 
surfaces of the telescope. Finally, the knowledge of the dust 
concentration at different sizes is also important for the proper 
design of the sealing of the telescopes and instruments moving 
mechanical parts. 

The first detailed analysis of more than 10 years of meteoro- 
logical data obtained using the TNG meteorological station at 
ORM can be found in the two recent papers, Lombardi et al. 
(120061 1 (hereafter Paper I) and Lombardi et al. (120071 1 (hereafter 
Paper II). Paper I shows a complete analysis of the vertical tem- 
perature gradients and their correlation with the astronomical 
seeing. Instead, Paper II shows an analysis of the correlation 
between wind and astronomical parameters, as well as the 
overall long-term weather conditions at ORM. Differences in 
the ORM microclimate demonstrated in a detailed comparison 
between synoptical parameters taken at three different locations 
in the observatory on a 1000 m spatial scale. Furthermore, 
ORM was shown to almost be dominated by high pressure and 
characterised by an averaged relative humidity lower than 50%. 
Extended dust clouds periodically appear above the ORM 
complex. Murdin (1985) shows that dust collected at the ground 
level at ORM is composed of small round quartz aggregates 
having diameters ranging between 10 to 60 fim diameter. The 
optical properties of this dust have been investigated by several 
authors (Stickland [19871 Murdin [19851 Jimenez et al. [19981 
all founding a flat wavelength dependence of the extinction; in 
particular, Murdin points out that Saharian dust's size is much 
bigger than the standard photometric V band wavelength. 
Big plumes of dust during the Calima days are also visible 
using satellites probes. In Siher et al. (2004) an analysis of 
5200 data covering the years 1978-1993 shows a correlation 
between astronomical extinction and the Total Ozone Mapping 
Spectrometer (TOMS) satellite derived extinction as well as a 
seasonal trend above the Canaries with a significant increase 
in dusty days at the end of July. Siher et al. ( 120041 ) also shows 
an average of 55 dusty days from the beginning of June to the 
beginning of September (about half of the summer nights). In 
this paper only the aerosol data corresponding to a significant 
absorption were used. The data since 1978 seems to indicate 
that the dusty period is drifting during the years. An indepen- 
dent analysis, based on Carlsberg Automatic Meridian Circle 
(CAMC) telescope data, has been published by Guerrero et al. 
([19981 . 

The present paper has the goal of improving the knowledge 
of the aerosol-particle-size distribution and the effects on astro- 
nomical observations in La Palma. In particular, we have anal- 



ysed for the first time the influence of the dust collected at the 
level of TNG primary mirror (Ml) as a function of its size to 
the aim: a) to understand if local measurements may be taken 
as good indicators of the upper conditions in the atmosphere (a 
dust ground condition has reproduced a similar situation in the 
atmosphere?); b) to give an accurate extinction study selecting 
among different sizes of dust; c) to use this information as a pri- 
mary tool for a fast warning of the unclean environment. 
The paper is organised as follows. Section [2] describes of the 
adopted particles counter. Section [3] describes the adopted pro- 
cedures to analyse the collected data. Section [4] defines and re- 
ports number of storm events. A comparison between local ex- 
tinction and total visual extinction obtained from the database of 
the CAMC is discussed in Section[5] 

Using TOMS satellite archive, we have extracted the Aerosol 
Index from TOMS/EarthProbe, which provides independent, 
semi-quantitative indications of the dust content in the high at- 
mosphere. Discussions of aerosol and TOMS data are in Section 
[6] conclusions can be found in Section[7] 

2. The dust monitor 

The selection of air-particles counting equipment is primary trig- 
gered by the particle sizes to be monitored to obtain their esti- 
mated amount. The parameters used for this selection were the 
sizes of the dust, the flow rate, and the background noise counts. 
Each measured size of the dust corresponds to the number of 
available channels of the counter. The sensitivity of the detector 
is given by the smallest size of particle that the sensor is able 
to detect. The flow rate is the amount of the flux that the sensor 
can receive through the sample volume, while the background 
noise is the instrument dark current value. We have addressed 
our choice to an easy-to-use instrument, hand-held, reliable, and 
very sensitive. Moreover, a particle counter based on the opti- 
cal detection technique that can count particles above 5.0 //m 
diameter was selected. The reason for such sensitivity is that 
we expect the number of dust particles to be less than the to- 
tal number of atmospheric particles and we expect them to occur 
mostly on scales near 1 /mi with a typical concentration range of 
a few tens per cubic centimetre. We therefore chose the particle 
counter Abacus TM301 (made by Particle Measuring System, 
Inc.), which uses a laser scattering technique for environmental 
ambient air analysis. Abacus TM301 is a compact and portable 
system designed to measure the purity of close environments, 
like a clean room, the counter is a system to measures the den- 
sity of particles in open air. This system is based on a filter 
through which the air is forced to flow. Different kinds of fil- 
ters are used sequentially to measure different sizes of particles 
centered at 0.3, 0.5, 1.0, and 5.0 yum according to the estimated 
dimension of dust and sand of the Sahara desert as described 
by Murdin dl9851 l. A small pump placed inside the dome sucks 
the air and, through scattering the light from a laser diode, mea- 
sures the number of particles. The four sizes of particles can be 
measured in both integrated and differential modes. Although the 
counter is mounted inside the TNG dome and at the level of tele- 
scope, aim long silicon pipe through the wall feeds the pump 
to external air at 13 m above the ground. This is a needed pre- 
caution because, even if we intend to use the counter to monitor 
the external air, it must be kept in a closed site, safe from ad- 
verse weather conditions. Unfortunately, we cannot distinguish 
between the type of particles, but only their size. For instance, 
it is not possible to distinguish between water vapour and dust. 
For this reason it is necessary to stop the monitoring if relative 
humidity increases to the condensation point (typically > 85%). 
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Fig. 1. Distribution of the atmospheric particles as a function of the different size, as measured by the Abacus TM301 during the 
years 2001-2006 (1-2 weeks each month depending on weather conditions). The gaps correspond to interruptions in the counter 
operations or to rejected values due to relative humidity > 85%. 



Table 1. Main characteristics of the Abacus TM301 dust moni- 
tor. 



size sensitivity 
input flow rate 
size channels 
light source 
concentration limits 
zero counts 
output 
data storage 
sample time 
purge time 
dalay time 



0.3 yum at 50% counting efficiency 

0. 1 Cubic Foot Minute [CFM] 

0.3, 0.5, 1.0, 5.0 fim 

laser diode (A = 780 nm) 

> 10 6 per dm instantaneous 

< 1 per CF 

RS-232 

500 sample, rotating buffer 
selectable from 1 sec to 99 minutes 
selectable from 1 sec to 99 minutes 
selectable from 1 sec to 99 minutes 



Table[T]summarises the instrumental basic performances already 
reported in Porceddu et al. (2002) and also reported here for 
completeness. 



3. Data analysis 

The dust monitor measures outdoor dust at the same level as 
to TNG Ml to have a monitoring as close as possible the Ml 
local conditions. Data coming from Abacus TM301 are consid- 
ered not as instantaneous values but as mean values after an in- 
tegration time of 10 minutes. Before each measurement, a 5- 
minute long reset is used to clean the sensor. Abacus TM301 
stores about 500 sets of data in memory, and when the mem- 
ory is full, it uses an RS232 serial common port to transfer data 
to the PC for further analysis. Up to now the database collects 
measurements from August 2001 to December 2006, and it is 
used consecutively for 1-2 weeks each month depending to the 
weather conditions. In particular the counting is stopped if rela- 
tive humidity > 85%. In Figure[T]we report the full dust-count 
distribution for each particle size as measured during the years. 
The gaps in the figure corresponding to zero counts are caused 
by interruptions in the counter operations or by rejected values 
due to relative humidity > 85%. 

To check the presence of a seasonal variation, we split the dust 
database in two epochs, defining wintertime the months from 
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Fig. 2. Seasonal distribution of the monthly dust at ORM: (top) 
0.3 pm dust; (middle - top) 0.5 pm; (middle - bottom) 1.0 fim 
dust; (bottom) 5.0 /mi dust. Different years are indicated by dif- 
ferent symbols. 

October until March and summertime the months from April un- 
til September. In this analysis we define the dust storm event by 
each dust count having a value a few cr-levels over the monthly 
mean values. The background is evaluated using the cr-clipping 
algorithm also described in Huber ( 1981 ) and Patat ( 120031) . 
The particle counts (N,) were assumed to have a Poissonian dis- 
tribution. For each month we computed the median value of the 
counts (MED), then for each different size of dust we calculated 
the median absolute deviation (MAD) defined as the median of 
the distribution 

\Ni-MED\. (1) 

We set the cr parameter to the value cr = 1 .48 MAD (accord- 
ing to Huber 1981 ), where 1.48 is the ratio between the standard 
deviation and the MAD under the assumption of a Gaussian dis- 
tribution. During the iteration we rejected counts having 

\Ni - MED\ > Kcr. (2) 

The procedure for estimating the dust background was iterated 
twice. In the first iteration we rejected the points exceeding +3<x 
with respect to the computed monthly median. In the second 
iteration, we rejected points exceeding 2cr and -3cr so that the 
background distribution is uncontaminated by any other dust 
peaks corresponding to dust storms. 
We have found that, in connection with dust storms, counts 

Table 2. Dust background content in [N rrT 3 ] at ORM in winter- 
time, summertime, and in the entire annual cycle. 





wintertime 


summertime 


annual 


0.3 yum 


1.27 ■ 10" 


4.38 • 10" 


2.98 • 10" 


0.5 yum 


1.21 ■ 10 s 


3.76 • 10 s 


2.53 ■ 10 5 


1.0 yum 


0.51 ■ 10 s 


1.51 • 10 5 


1.02 • 10 5 


5.0 yum 


0.66 • 10 3 


1.47 ■ 10 3 


1.09 - 10 3 



Fig. 3. Dust storm event of 2002 December 25 and 26: 0.3 pm 
(solid), 0.5 //m (dots), 1.0 pm (short-dashes) and 5.0 /im (long- 
dashes). 

of particles having diameters of 0.3, 0.5, and 1.0 pm increase 
by 2 orders of magnitude with respect to the background. 
Unfortunately 5.0 pm particles often present a typical fluctua- 
tion between 1 and 2 orders of magnitude, so it is difficult to 
distinguish dust storms from the typical dust background. For 
this reason an increase of 3 orders of magnitudes does allow us 
to distinguish the 5.0 pm peaks from the background. 
Figure [2] shows the median monthly distribution of the particles 
as a function of the month. It is evident that the level of back- 
ground depends on size. The distribution shows the increase in 
the counts during February-April and July-September, while 
May and June are the months with the lower level of dust. 
Table [2] shows the median dust content as computed in the two 
seasons and in the entire annual cycle. The table clearly shows 
that small particles are dominant both in wintertime and in sum- 
mertime. It is not surprising because heaviest particles decay 
quickly, while lighter particles flow in the atmosphere for longer. 



4. Dust storms 

As mentioned, dust storm is the event where the content of dust, 
given in counts per cubic metre [N trr 3 ], increases with respect 
to the typical dust background content in clear days by 2 orders 
of magnitude or more for 0.3, 0.5, and 1.0 /im dust, and 3 orders 
of magnitude or more for 5.0 pm dust. Figure [3] shows a more 
detailed situation of a dust storm that occurred in December 
2002. This event persisted for two days, on 2002 December 25 
and 26, and it shows an increase of several orders of magnitude 
in the counts. 

In a deeper analysis we noticed that the dust storms are 
typically 3-4 days long, but sometimes it is possible to have 
some short storms 2 days long or, more rarely, long-lasting 
storms persisting 5-6 days. We correlated the wind direction 
with the dust-storm counts. Figure [5] shows the wind roses 
obtained during the storm events as a function of particle size. 
The wind direction is recorded by the TNG meteorological 
tower and data were already analysed in Paper II. There is a 
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clear correlation between dust storms and wind direction. Only 
the winds blowing from south-west, i.e. from the edge of the 
Caldera de Taburiente, carry a large amount of dust particles. 
We confirm that the best photometric conditions occur when the 
wind blows from the north. 

In Table [3] the total number of events for both wintertime and 
summertime in the period from August 2001 to December 2006 
is reported. Particles having 1.0 and 0.5 fim size show the larger 
number of events. Table [3] shows that dust storms occur with 
the same statistic in both wintertime and summertime. Whittet 
et al. ( 119871 ) demonstrated that 95% of the aeolian deposits 
collected at La Palma have diameters between 1 and 100 pm. 
To better understand this analysis we calculated the median 
value of the dust content in case of storms for each measured 
size, and to include more massive grains we extrapolated the 
counts to 10.0 fim. Following Whittet et al. (119871) . we assume 
a power-law distribution r 2 N(r) for particles > 1 .0 /mi. Because 
we cannot distinguish between seasons, we decided to consider 
this background as a constant in both wintertime and summer- 
time. Under these assumptions we found that 10.0 //m particles 
have a typical background of ~10 N irT 3 which is a lower 
concentration than the other particles we measured. The same 
power law permits us to extrapolate the concentration to the 
10.0 fim size in case of dust storms. We obtained a concentration 
of ~4- 10 4 N m~ 3 . Results of this analysis are plotted in Figure[4] 



5. Aerosol atmospheric extinction in B, V, and / 

The optical properties of the dust in La Palma have been investi- 
gated by several authors who found that aerosol extinction does 
not have a clear dependence on the adopted wavelength. In par- 
ticular, Stickland (1987) pointed out that extinction during dust 
storms is relatively stable on a time scale of 15 minutes. No anal- 
ysis exists for the effects of the different dust sizes on astronom- 
ical extinction. 

In this paper we want to compare the differential extinction com- 
puted using our measured dust with simultaneous nighttime as- 
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Fig. 4. Distribution of the median dust counts during dust storms. 
Particles >1.0 yum are treated as if they follow an r 2 N(r) power 
law. 







Fig. 5. Wind roses at TNG during dust storm events. 



tronomical extinction obtained from the CAMC database, which 
is a very useful tool for quantitatively studying dust effects. We 
evaluate aerosol atmospheric extinction (k) for each size of dust 
at 3 different wavelengths: 435, 550, and 780 nm correspond- 
ing to the central wavelengths of the standard photometric bands 
B, V, and /. We believe that using the differential dust counter 
may help us evaluate the different extinction contributions from 
the different dimensions of the grains as a function of the wave- 
length. 

The aerosol extinction is obtained by applying the model based 
on the Mie theory and has also been used by Mathis et al. 
(fT977l> Whi ttet et al. JT9871 ), Jimenez et al. ( fT998t , and Marley 
et al. d 19991 ). The Mie theory assumes the particles are homo- 
geneous spheres of radius r. Following Patterson ( 119771 ) defini- 
tions, the scattering coefficients for a Mie solution are linked to 
the size and composition of the particles through the parameter 
x = 2nr/A. For an incident light at wavelength A, the extinction 
coefficient k due to particles of radius r is defined by the formula 



kxs = Qext(n, r, A) nr 2 N, 



(3) 



where the term Q ex t(n, r, A) is the extinction efficiency factor 
computed for Mie's particles, while N is the number of parti- 
cles of radius r in a cubic centimetre. With this definition, kx, r is 
expressed in magnitudes per centimetre [mag cirT 1 ] and is corre- 
lated with the local density of the dust content in the atmosphere 
above the TNG dome. 

In the equation, the extinction efficiency factor Q ext is a com- 
plex function depending on the particles radius r, the absorbed 

Table 3. Number of storm events at ORM in wintertime, sum- 
mertime, and in the entire annual cycle in the period from August 
2001 to December 2006. 





wintertime 


summertime 


total 


0.3 yum 


10 


8 


18 


0.5 yum 


15 


12 


27 


1.0 /im 


21 


22 


43 


5.0 yum 


9 


7 


16 
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Fig. 6. Extinction efficiency factor as a function of the particles 
radius in B (dots), V (solid), and / (short-dashes). 

wavelength A and the complex refractive index n. In the defini- 
tion (Patterson 1981 ), the real part of n is the ratio between light 
speed in a vacuum and the light speed in particles, while the 
imaginary part is related to the Bouguer absorption coefficient 
k through the formula n !m = kA/Aji. Following Patterson (119771 1, 
Carlson & Benjamin ( 1980), Fouquart et al. ( 1987), and Jimenez 
et al. (fT998l , we assigned n = 1.55 - 0.005/ (i = V^I) as the 
typical refractive index for Saharan dust aerosols. 

We computed the extinction efficiency factor Q ext for the 
adopted photometric bands using the Mie scattering code of 
Wiscombe ( 119801 ). Figure [6] shows the plot of Q ext as a func- 
tion of the particle radius for B, V, and / bands, while in Table [4] 
we report the computed values of Q ext for 0.3, 0.5, 1.0, and 5.0 
fim dust in the three bands. We show that for r > 1 .0 fim, Q ext is 
almost constant in the three bands, confirming that for massive 
grains the extinction is colour-neutral as pointed out by Jimenez 
et al. (TT9981 . 

Applying Eq. (0, we obtained the local aerosol atmospheric ex- 
tinction in [mag cirT 1 ] caused by different grain sizes in typical 
dust background conditions (^med) and in the presence of dust 
storm events (k^ usty ). 

Table [5] reports the Icmed in B, V, and / in wintertime, sum- 
mertime, and in the entire annual cycle using data as reported 
in Tableland for the extrapolated value of dust background at 
10.0 fim. As shown in the table, on clear days 0.3 fim dust gives 
the highest contribution to the global extinction, in particular in 
summertime. The extinction due to 0.3 /mi dust is lower at higher 
wavelengths, while for 0.5 fim dust it becomes stronger at higher 

Table 4. Values of the extinction efficiency factor Q ext . 





B 


V 


/ 


0.3 fim 


4.148 


4.334 


2.783 


0.5 fim 


1.946 


2.756 


4.137 


1.0 fim 


2.367 


2.355 


2.205 


5.0 fim 


2.142 


2.135 


2.140 


10.0 fim 


2.072 


2.084 


2.099 



Table 5. Seasonal and annual local aerosol atmospheric extinc- 
tion computed in dusty background conditions^. 







wintertime 






kMED(B) 


k M ED(V) 


klAED (I) 


0.3 fim 


1.49 


1.56 


0.99 


0.5 fim 


0.19 


0.26 


0.39 


1.0 fim 


0.38 


0.38 


0.35 


j.V fim 


n i 1 
U. 1 1 


n 1 i 


n 1 1 
U.l 1 


inn, im 
iu.u fxm 


< U.U1 


< u.ut 


< u.ut 






summertime 






kMED(B) 


k M ED(V) 


klAED U) 


0.3 fim 


5.14 


5.37 


3.45 


0.5 fim 


0.58 


0.81 


1.22 


1.0 fim 


1.12 


1.12 


1.05 


5.0 fim 


0.25 


0.25 


0.25 


10.0 fim 


<0.01 


< 0.01 


<0.01 






annual 






kMED(B) 


k M ED(V) 


khAED (I) 


0.3 fim 


3.49 


3.65 


2.34 


0.5 fim 


0.39 


0.55 


0.82 


1.0 fim 


0.76 


0.76 


0.71 


5.0 fim 


0.18 


0.18 


0.18 


10.0 fim 


<0.01 


< 0.01 


<0.01 



Values are given in [10 8 mag cm ']. 



wavelengths, in particular in the summertime. Finally, the effect 
of the grains starting from 5.0 yum is almost constant in the three 
bands. We found that the contribution of the extinction due to 
10.0 fim particles is negligible. 

To understand the aerosol contribution on the total astronomical 
extinction, we need to evaluate the distribution in altitude of the 
particles. We compare our local extinctions in V to simultaneous 
astronomical extinctions in V from the CAMC database (kcAMc) 
under the assumption that kcAMc is mainly due to our measured 
dust. We assume 

kaer(X) = ^,0.3 + ^,0.5 + ^,1.0 + ^i,5.0 + ^.t, 10.0 (4) 

Because the CAMC-extinction database is only in V band, we 
compare our V band aerosol extinctions. For each k aer (V) we 
assume that particles are distributed in the atmosphere at an alti- 
tude h, in order to have 

kaeriV) h = kcAMc(V). (5) 

We assume again that dust is uniformly distributed in altitude in 
a column of atmosphere having a base of 1 cm 2 and height h 
as obtained from Eq. @. Figure [7] shows the distribution of the 
altitudes above the ORM obtained from each computed mean 
dust count. The maximum of the histogram shows the altitude 
reached by the dust. We found that the peak of the altitude dis- 
tribution of the dust is about 2.5 km above the telescope. This 
value agrees very closely with studies based on independent ra- 
diosonde data (Hsu 1999) and indicates that the settlement pro- 
cess of the dust from the Sahara to La Palma, in the size range 
0.3-10.0 fim, is not relevant, in agreement with the models of 
Murphy (1990). We conclude that, under the assumption of a 
uniform distribution of the dust, the particles are distributed in 
the first 2.5 km, on average. We then recompute the aerosol at- 
mospheric extinction using Eq. (0, integrating on a column of 
atmosphere of 2.5 km. We assumed this obtained altitude to be 
equivalent to 1 airmass and the final extinction is therefore given 
in [mag airmass -1 ]. 
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Fig. 7. Distribution of the dust altitudes from ORM altitude for 
Eq. ©. 



The kx, r for each particle size in typical dust background con- 
ditions and in typical dust storm conditions are now recomputed 
by integrating on a 2.5 km column of atmosphere above the TNG 
mirror. Figure|8]shows the aerosol atmospheric extinction k^ r in 
[mag airmass -1 ] and integrated on a column of atmosphere of 
2.5 km in the entire annual cycle as a function of the size of the 
dust (cfr. Tableland Figure|4| for clear and dusty days. 
We show that on clear days, the extinction is dominated by 0.3 
fj.m particles. Table [6] reports the percentage of contribution by 
different size of dust in each filter on both clear and dusty days. 
As shown in the table, on clear days 0.3 /mi particles are respon- 
sible for more than 70% of the total aerosol extinction in B and 
V, while their contribution in I is about 60%. The biggest parti- 
cles (>1.0 /mi) dominate in dusty days. It is interesting to note 
the different behaviors of I band during dust storms with respect 
to the B and V bands. Bigger particles show a contribution to the 
total extinction of more than 70% in all the bands (see Table[6]l. 
Following Eq. (@}, the expected total local aerosol atmospheric 
extinction has to be calculated as the sum of the contributions 
of each particles size. On dusty days this extinction is typically 
~0.2 mag airmass -1 (see Figure [8]l and can rise to ~1.0 mag 
airmass -1 during very strong dust storms (see Figure |5). Figure 
[9] shows k aer (V) calculated on dusty days compared to simulta- 
neous astronomical extinctions in V band as extracted from the 
CAMC database. Figure [9] shows a good correlation, the dotted 
line is the linear fit at a confidence level of 0.8 and having a slope 
of 0.77. We can conclude that it is possible to have an estimation 
of astronomical extinction, starting from a measure of dust count 
and by an appropriate calibration. 

To evaluate the increase in the local aerosol extinction in con- 
nection with dust storm events, we also evaluated the differences 
between (i) the aerosol extinctions calculated on each day of a 
dust storm and (ii) the extinction calculated in the typical dust 
background conditions of each defined seasonal epoch: 

Ak(A) = k dusly (X) - k MED (X). (6) 

Figure [l0]plots the mean increases calculated for each dust size 
in B, V, and /. In the figure we have also taken the extrapolated 



Fig. 8. Median aerosol atmospheric extinction in B, V, and / 
in typical dust background conditions and in typical dust-storm 
conditions for each particle size. The biggest particles are domi- 
nant with respect to sub-micron particles. 

extinction due to 10.0 //m particles into account. The typical in- 
creases are between 0.015 and 0.065 mag airmass . 

6. Summertime aerosol atmospheric extinction and 
TOMS aerosol index 

Since from the extinction data alone it is not possible to dis- 
tinguish cirrus clouds from dust clouds, we have combined both 
satellite and ground-based data comparing the daily k aer (A) com- 
puted in summertime in the three photometric bands and the tro- 
pospheric TOMS aerosol index (AI) that characterises the lo- 
cal aerosol concentration using satellite-based reflectivity mea- 
surements. The TOMS data are set with the data closest to the 
ORM within a 1 x 1 deg latitude and longitude box. TOMS AI 
is selected on the same days as k aer {X). Following Siher et al. 
(120041 ) we have set our threshold at AI > 0.7. Because TOMS AI 
measurements above the Canaries are done in daytime, we ob- 
tained the nighttime AI by interpolating the values of two con- 
secutive days. Satellite AI is obtained by the reflectivity effect. 
Because both dust and clouds may show reflectivity even at dif- 
ferent threshold, we rejected AI points when sky reflectivity > 

Table 6. Relative contribution of each particle size as percentage 
of the total aerosol atmospheric extinction k aer . 
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Fig. 9. k aer versus kcAMc in the case of dust storms. The linear fit 
has a confidence level of 0.8. 

15% (Bertolin 2005 ) to be sure to have only the dust effect. 
Once more we split the data between wintertime and summer- 
time, and we noticed the absence of AI points in wintertime in 
connection with our local aerosol extinction calculations. Most 
of those points are the ones previously rejected for high sky re- 
flectivity. In fact, in wintertime the sky is often contaminated 
by high reflectivity caused by thin cirrus or snow at the ground. 
Furthermore, in wintertime the dusts are often at low altitudes, 
too deep in the atmosphere to be revealed by the satellite. 
FigureQT]shows the plot of the local aerosol atmospheric extinc- 
tion in summertime versus the AI in V band (B and / show sim- 
ilar behaviours). In the plot we can distinguish three areas. The 
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Fig. 10. Typical local aerosol extinction increases Ak(A) in [mag 
airmass -1 ] in B (top), V(center) and / (bottom) in connection 
with dust storm events at ORM. The filled circles represent the 
expected Ak(A) due to 10 yum particles. 



Fig. 11. Summertime local aerosol atmospheric extinction ver- 
sus TOMS Aerosol Index in V band. The correlation calculated 
with the Spearman's test between the points inside the zone de- 
limited by the solid and the short-dashed lines has confidence 
level > 0.8. 

locus between the solid and short-dashed lines include all points 
having a correlation between k aer (A) and AI with a confidence 
level > 0.8 (Spearman's test). The points below the short-dashed 
line show high AI but low local extinction. This high reflectivity 
is caused by the Saharan dust suspended at a very high level that 
is not detectable at the TNG mirror level. The (very few) points 
above the solid line show high local extinction and low AI. We 
ascribe the cause to local, recycled dust that remain suspended 
at low altitude and cannot be revealed by the satellite. 
To support this affirmation, we computed the two wind roses 
for the data subsets relative to those days on which we believe 
that the extinction is due to local recycled dust, and on those 
days when the high reflectivity comes from the Saharan dust 
suspended in the high atmosphere. We found that, in the case 
of local recycled dust, the winds come prevalently from W-SW, 
while in the case of dust suspended in the high atmosphere the 
prevailing direction is from S-SW. 

We can conclude that observations may often be affected also 
by local recycled dust and not only by Calima events. More data 
may help us this point. Taking all the reasons mentioned above 
into account we conclude that in site testing, both local dust 
monitors and data from satellites have to be used together for 
a complete investigation of the effects of the aerosols. 

7. Conclusions 

The high-altitude transport events of dust from desert areas 
have been recognised the major atmospheric events on the 
Earth. The resulting pollution has astronomical observations 
and instrumentation (i.e. the main mirror reflectivity) effected 
even at very great distances, in particular at the Observatorio del 
Roque de Los Muchachos, on La Palma Island. 
In this paper we have presented an analysis of the properties of 
the local atmospheric dust measured at the height of the TNG 
primary mirror, the correlations of the dust peak events with the 
wind direction, and the total atmospheric extinction using the 
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5-year long data archive of the TNG dust monitor. This particle 
counting system is able to detect particles of 0.3, 0.5, 1.0, and 
5.0 /im, thereby according with the estimated dimension of 
dust and sand of the Sahara desert. Typical annual average 
background values are 3 • 10 6 , 2.5 • 10 5 , 10 5 , 10 3 particles 
per cubic metre, respectively, in the four channels. During 
dust storm events, the particle concentration increases to three 
orders of magnitude. Using a power-law distribution we also 
extrapolated the contents of 10.0 //m dust on dusty days. We 
estimated an amount of 4 • 10 4 particles per cubic metre for 
massive grain. 

After considering two seasonal epochs, we found the number 
of dusty events to be comparable, but they reach a somewhat 
higher concentration in summertime, when the background is 
also higher by a factor 3-4. In both seasons the concentration 
of small dust particles (0.3 yum) is dominant with respect to 
the other size dust particles. A study of dust-storm events 
versus wind direction has shown that winds coming from the 
southwest, i.e. from the edge of the Caldera de Taburiente, 
typically carry a large amount of dust particles. Dust storms 
are 3-4 days long, but sometimes it is possible to have some 
shortly storms 2 days longer or, more rarely, long-lasting storms 
persisting 5-6 days. 

A comparison between the storm duration and the TOMS 
Aerosol Index shows that short-storms typically are not linked 
with a high Aerosol Index. This suggests that short-storms can 
be generated by local reprocessing of currents carrying dust in 
suspension at low altitudes such is the TNG mirror level. We 
also noticed that dust storms occur not only in summertime 
(Saharan Calima events), but that several events occur during 
other seasons as effect of local currents recircle . 
In clear days the aerosol extinction is dominated by 0.3 fim 
particles (~70%), while in dusty days the strongest contribution 
is from bigger particles having size >1.0 fim (>70%). On dusty 
days the extinction is typically ~0.2 mag airmass -1 , and it can 
rise to ~1.0 mag airmass -1 in the presence of very strong dust 
storms. 

In dusty days dust V extinction shows a positive correlation 
with simultaneous astronomical extinctions values from CAMC 
database (c.l. 0.8). This suggests that the dust detected at TNG 
affects the astronomical observations in a very wide area. 
Assuming a uniform distribution of the dust in the observing 
column, we found that the dust extends about 2.5 km above the 
telescope, corresponding to an height of 5 km above the see 
level. This confirms that the settlement process of the dust in the 
size range 0.3-10.0 fim, is not relevant. 



Huber, P. J. 1981, Robust Statistics, ed. J. Wiley, & Sons, New York 

Hsu, N. C, Herman, J. R., Torres, O., Holben, B. N., Tanre, D., et al., 1999, J. 

Geophys. Res., 104, 6269 
Lombardi, G., Zitelli, V., Ortolani, S., & Pedani, M. 2006, PASP, 118, 1198 

(Paper I) 

Lombardi, G., Zitelli, V., Ortolani, S., & Pedani, M. 2007, PASP, 119, 292 (Paper 
II) 

Marley, M. S., Gelino, C., Stephens, D., Lunine, J. I., & Freedman, R. 1999, ApJ, 
513,879 

Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Murdin, P. 1985, Vistas Astron., 28, 449 

Murphy, J. R., Toon, O. B., Haberle, R. M., & Pollack, J. B., 1990, J. 

Geophys.Res., 95, 14629 
Patat, F. 2003, A&A, 401, 797 
Patterson, E. M. 1977, Appl. Opt., 16, 2414 
Patterson, E. M. 1981, Science, 211, 836 

Porceddu, I., Zitelli, V., Buffa, F, & Ghedina, A. 2002, SPIE, 4844, 358 
Siher, E. A., Ortolani, S., Sarazin, M., & Benkhaldoun, Z. 2004, SPIE, 5489, 138 
Stickland, D. J., Lloyd, C, Pike., C. D., & Walker, E. N. 1987, The Observatory, 
107, 74 

Wiscombe, W. J. 1980, Appl. Opt., 19, 1505 

Whittet, D. C. B., Bode, M. F, & Murdin, P. 1987, Vistas Astron., 30, 135 



Acknowledgements. The authors acknowledge the anonymous reviewer for the 
very useful comments and the TNG staff for the availability of the data. The 
authors thank also Guido Barbara of Department of Astronomy of University 
of Padova for helpful comments on Mie's theory and Chiara Bertolin of 
Atmospheric and Climate Science Institute of National Research Council of 
Padova for the help in using TOMS/EarthProbe data. 



References 

Bertolin, C. 2005, private communication 

Carlson, T. N., & Benjamin, S. G. 1980, J. Atmos. Sci., 37, 193 

Fouquart, Y., Bonnel, B., Chaoui Roquai, M., Santer, R., & Cerf, A. 1987, J. 

Climate Appl. Meteor., 26, 28 
Ghedina, A., Pedani, M., Guerra, J. C, Zitelli, V., & Porceddu, I. 2004, SPIE, 

5489, 227 

Guerrero, M. A., Garci'a-Lopez, R. J., Corradi, R. L. M., Jimenez, A., Fuensalida, 

J. J., et al., 1998, New Astron. Rev., 42, 529 
Jimenez, A., Gonzalez Jorge, H, & Rabello-Soares, M. C. 1998, A&AS, 129, 

413 



